Evolution of the nuclear and magnetic structures of TlFei 6Se2 with temperature 
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The evolution of the nuclear and magnetic structures of TlFei.6Se2 was determined in the tem- 
perature range of 5-450 K using single crystal neutron diffraction. The Fe layers in these materials 
develop a corrugation in the magnetically ordered state. A canting away from the block checkerboard 
magnetic structure is observed in the narrow temperature range between approximately 100 and 
150 K. In this same temperature range, an increase in the corrugation of the Fe layers is observed. 
At lower temperatures, the block checkerboard magnetic structure is recovered with a suppressed 
magnetic moment and abrupt changes in the lattice parameters. Microstructure analysis at 300 K 
using atomic-resolution Z-contrast scanning transmission electron microscopy reveals regions with 
ordered and disordered Fe vacancies, and the iron content is found to be uniform across the crystal. 
These findings highlight the differences between the alkali- metal and thallium materials, and indi- 
cate competition between magnetic ground states and a strong coupling of magnetism to the lattice 
in TlFei.6Se2. 

PACS numbers: 74.70.Xa, 75.25.-j 



The recent discovery of superconductivity at tem- 
peratures greater than 30 K in AxFey^e2 is generat- 
ing much effort to understand the relationship between 
superconductivity, magnetism, and structure in these 
materials.^ These compounds crystallize in a vari- 
ant of the ThCr2Si2-type structure, with vacancy or- 
dering in the iron layer yielding a y/ba x ^/5a super- 
cell. Although TlFei.6Se2 is a magnetic insulator at low 
temperatures,^^ this compound is chemically and struc- 
turally similar to recently discovered superconducting 
compositions where part or all of the Tl is replaced by an 
alkali metal such as K, Cs, or Rb.^ ^ For instance, simi- 
lar to the superconducting compositions, vacancy order- 
ing occurs below Ty ^ 460 K in TlFei.6Se2, and a block 
checkerboard antiferromagnetic magnetic order emerges 
below Tat - 430 K.l^ There are significant differences, 
however, between the superconducting compositions and 
TlFei.6Se2. 

For TlFei.6Se2 crystals grown from the melt, the Tl 
site is fully occupied, and the total Fe content shows lit- 
tle variation regardless of the starting composition. It is 
likely, however, that the degree of vacancy order depends 
on the starting composition and the processing condi- 
tions. The superconducting compositions appear to be 
more flexible, and the best superconducting properties 
are obtained when the alkali metal and iron site occu- 
pancies are both near 80%, e.g., Ko.8Fei.6Se2.™™^ 
One distinct advantage of the TlFei.6Se2 crystals is that 
they are stable in air for many hours and do not exhibit 
any tendency toward phase separation into Tl-rich and 
Tl-poor regions, which appears to occur when Tl is re- 
placed by an alkali metal.^^^^ 

Our previous work on TlFei.6Se2 revealed two unex- 
pected phase transitions, at Ti ^ 140 K and T2 ^ 100 K, 
which were evident in resistivity, magnetic susceptibility 
and heat capacity data.^^ Upon cooling, the magnetic 



order parameter was observed to reach a maximum at 
~ Tl , below which it decreased with a particularly sharp 
decrease just below T2. These transitions did not cor- 
respond to a change in the crystal symmetry or vacancy 
ordering, and understanding their physical origin was the 
motivation for this detailed structural study. 

In the present work, the nuclear and magnetic struc- 
tures in the vicinity of the two phase transitions Ti and 
T2 are examined in detail using single crystal neutron 
and X-ray diffraction. The possibility of phase sepa- 
ration into fully-ordered and disordered regions is ex- 
amined by atomic-resolution Z-contrast scanning trans- 
mission electron microscopy (Z-STEM). Such phase sep- 
aration has been re ported in crystals and films of 
K,Fe2-,Se2,>niIlIinin]and in Rbo.gFei.eSeaP 

Single crystals of TlFei.6Se2 were grown from a melt 
of nominal composition TlFe2Se2 using a modified Bridg- 
man technique as discussed in Ref. [161 The crystals were 
observed to be homogeneous by energy dispersive spec- 
troscopy, with a resolution on the order of 50 jiia. Sin- 
gle crystal neutron diffraction measurements from 5 K to 
450 K were made at HB-3A four-circle diffractometer at 
the High Flux Isotope Reactor at the Oak Ridge National 
Laboratory. Complete neutron diffraction data sets were 
collected at 5K, 90 K, 115 K, 140 K, 175 K, 250 K, 340 K, 
and 450 K. The neutron wavelength of 1.536 A was 
used with a bent perfect Si-220 monochromator.^ Sin- 
gle crystal X-ray diffraction was measured at 135 K and 
200 K with a Bruker SMART APEX CCD diffractome- 
ter with Mo-Ka radiation. The neutron diffraction data 
refinements are based on '^250 reflections and were com- 
pleted in the program Fullprof.l^ The data were refined 
using the A/5ax^/5a superlattice notation, and both nu- 
clear and magnetic structures were resolved. Refinement 
results from neutron data at selected temperatures are 
shown in Table HJ See Ref. 1161 for additional information. 
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FIG. 1. (color online) (a) HAADF image with the electron beam parallel to the a -axis of the y/5a x \/5a supercell. The 
bright lines correspond to the heavier Tl planes, (b) and (c) are digital diffraction patterns from the regions outlined in (a) by 
blue (disordered Fe vacancies) and red (ordered Fe vacancies) squares, respectively, (d) Fourier filter obtained by selecting the 
superstructure reflections indicated by circles in (c), with black regions corresponding to disordered regions. 



The microstructure of the crystals was examined at 
300 K using atomic-resolution Z-contrast scanning trans- 
mission electron microscopy (Z-STEM). Aberration- 
corrected high angle annular dark field (HAADF) imag- 
ing was performed in a Nion UltraSTEM equipped with 
cold field emission gun, third generation C3 and C5 
aberration corrector, and operated at 100 kV. Samples 
were thinned to nearly 20 nm via dual focused ion beam 
(FIB) slicing followed by Ar ion milling at a voltage of 
0.5 kV. FIB was employed as an alternative to conven- 
tional grinding and polishing to avoid introducing exces- 
sive strain in the sample. Low voltage ion milling was 
employed to remove FIB induced amorphization of the 
sample surface. The image was treated using a probe 
deconvolution algorithm to improve contrast and reduce 
noise. ^ 

Previous refinement of the neutron diffraction data re- 
vealed the magnetic moment to be l.T/i^/Fe at room 
temperature in TlFei.6Se2f^ This is smaller than ob- 
served in the analogous superconducting compositions, 
where the moment is ^S/i^/Fe below ~250K."^ Ad- 
ditionally, multiple, unexplained low-temperature mag- 
netic phase transitions were observed in the magnetic 
susceptibility, electrical resistivity and heat capacity of 
TlFei.eSes at Ti and T2 .^Hl Electron microscopy was em- 
ployed to determine if these results could be explained by 
phase separation into regions of differing vacancy order, 
which would likely have different magnetic behaviors 

Examination of the microstructure in TlFei.6Se2 using 
cross section Z-STEM reveals regions of vacancy order 
and disorder, as observed in Figure [l] Figure [l^ is an 
HAADF image with the electron beam directed along 
the a axis of the ^/5a x ^/ba supercell. A close in- 



spection of Fig. [T^ shows the existence of two regions 
yielding slightly different contrast and lattice periodic- 
ity. While regions like the one highlighted by the blue 
square labeled (b) yield fast Fourier transforms (FFT) 
that are consistent with a ThCr2Si2-type structure (Fig. 
[l]3), regions like the one framed by the red square show 
superlattice reflections associated with periodic minima 
of intensity in the Fe plane (see Fig. [l]^). This allows the 
ordered and disordered regions to be visualized, as in Fig. 
[l]i, which is a filtered image obtained by selecting only 
the frequencies yielding superstructure reflections indi- 
cated by the circles in Fig. [T]^. Figure [T]i highlights the 
spatial extent of those regions producing superstructure 
reflections, clearly delineating the ordered (light) regions 
from the disordered (dark) regions. Using Fig. and 
similar images (spanning areas of 64 x 64 nm^), a vol- 
ume fraction in the range 40 - 48% is obtained for the 
ordered regions. This quantification assumes that the 
ordered (disordered) regions extend for the entire thick- 
ness of the TEM specimen (nearly 20 nm) in the direction 
perpendicular to the image plane. 

Figure |2^ is a higher magnification HAADF image of 
the same sample. The center and right hand corners of 
the image have highly ordered vacancies, while the left 
hand corners are disordered. The green and blue rect- 
angles highlight unit cells in the ordered and disordered 
regions, respectively. 

Intensity line profiles along the Fe planes labeled 1 and 
2 in Fig. |2^ are reported in Fig. ^p. These reveal the lo- 
cation of vacant Fe sites as positions of the deep minima 
in intensity. Curve 1 exhibits deep minima in correspon- 
dence with the highly vacant Fe2 sites of the model in 
Fig. |2]i, as expected for a well ordered region with the 
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FIG. 2. (color online) (a) High magnification HAADF image 
with the electron beam parallel to a. The brightest planes 
are Tl planes. The nearest neighbor planes are Se planes, and 
the dimmest planes sandwiched between two Se planes are Fe 
planes. The green (blue) rectangle frames one unit cell, (b) 
Plot of the HAADF intensity along the Fe planes labeled 1 
and 2 in (a), (c) Temperature dependent integrated intensity 
of the Bragg peak (2 0), which is related to vacancy order, 
(d) Projection of the Fe layer of TlFei.6Se2 along the b axis; 
spheres represent Fe where the Fe site occupancies are shown 
by partial brown coloring with Fe moments indicated by red 
and blue arrows to highlight the spin blocks. 



^/5a X ^/ba supercell. Curve 2 is not entirely consistent 
with the model, as expected for a disordered distribution 
of vacancies. Analysis of intensity line profiles for the 
Fe planes over distances of ~35 nm suggests that the or- 
dered and disordered regions have on average the same 
number of vacancies. Therefore, the ordered and disor- 
dered regions have the same composition and differ only 
by the arrangement of the Fe atoms (vacancies). 

In summary, the Z-STEM data suggest a distribution 
of regions of '^10nm-30nm in size along the basal plane 
and 2-4 unit cells along the c -axis having a high con- 
centration of ordered Fe vacancies. The Z-STEM also 
reveals these ordered/disordered regions are crystallo- 
graphically coherent, which is consistent with the obser- 
vation of sharp diffraction peaks. 

The presence of ordered and disordered regions does 
not entirely explain the observed physical properties. 
Specifically, the observation of two magnetic transitions 



TABLE I. Refined structural parameters for TlFei.6Se2 from 
single-crystal neutron diffraction data at three temperatures; 
space group 14/ m (No. 87). All non- listed atomic coordinates 
are zero except for Fe2, which is located at (0, |, ^); data were 
refined with U^^^^U^^'^ and [/f/^ = U^f. 





5K 


115K 


250 K 


a(A) 


8.641(4) 


8.665(4) 


8.685 (4) 


c(A) 


13.993(7) 


13.957(7) 


14.002(7) 



XT12 
yTl2 
XFel 

ypei 

ZFel 
ZFe2 
ZSel 
XSe2 
ySe2 
ZSe2 



0.1954(13) 
0.3957(11) 
0.0928(7) 
0.1970(7) 
0.2476(6) 
0.727(5) 
0.3651(12) 
0.1986(10) 
0.3923(8) 
0.3565(3) 



0.1964(15) 
0.3959(15) 
0.0933(7) 
0.1964(8) 
0.2452(9) 
0.717(5) 
0.3638(14) 
0.1980(11) 
0.3925(10) 
0.3560(4) 



0.1943(18) 
0.3966(14) 
0.0926(8) 
0.1967(8) 
0.2485(8) 
0.722(5) 
0.3648(15) 
0.1993(12) 
0.3921(9) 
0.3557(4) 



eg (A 

jFel 



0.0184(16) 
0.0126(18) 

0.04(2) 
0.0185(14) 



0.0295(17) 
0.0164(20) 

0.07(3) 
0.0244(15) 



0.038(2) 
0.020(2) 
0.05(3) 
0.0261(18) 



Ri 
wRF^ 



0.0538 
0.122 
7.7 



0.0522 
0.123 
9.11 



0.0589 
0.132 
8.04 



cannot be understood as one transition occurring in the 
ordered region and the other transition occurring in the 
disordered region. Such a conclusion is inconsistent with 
the diffraction results discussed below. However, an in- 
teraction between these regions, most likely via lattice 
strain, may be responsible for the observed behavior. 

Single crystal neutron and X-ray diffraction both re- 
veal a vacancy ordering transition around 460 K. The 
temperature dependent integrated neutron intensity of 
the nuclear Bragg peak (2 0) is plotted in Figj2]3. The 
(2 0) reflection is solely due to ordering of the Fe va- 
cancies and fully disappears at ~ 460 K by extrapolating 
the observed data points. Upon cooling to ~ 400 K, the 
vacancy order saturates. Therefore, the Z-STEM data at 
300 K are expected to be representative of the microstruc- 
tures at lower temperatures as well. 

All of the diffraction data are well described by a 
>/5a X y/ba superstructure with a highly Fe-occupied po- 
sition Fel (on average 93% occupied) and a highly-vacant 
position Fe2 (on average 27% occupied), consistent with 
our prior report. The volume percentage of ordered and 
disordered regions estimated from the HAADF images is 
consistent with these average occupancies. The refine- 
ment results are summarized in Table [l] for T = 5K, 
115 K, and 250 K. 

The vacancy ordered Fe layer is shown in Figj2]i, with 
Fe occupancies proportional to the filling of the brown 
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spheres. The Fe moments in the block checkerboard 
structure are shown as red and blue arrows to highlight 
the blocks of spins, which are ferromagnetically aligned 
within a block and antiferromagnetic between blocks. 
The corrugation {S) between the spin blocks in Fig. |2]i 
is temperature dependent and will be discussed in detail 
below; 5 is defined as the distance along c that separates 
the spin blocks. 
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FIG. 3. In (a) and (c) the lattice parameters are shown to re- 
spond strongly to the magnetic transitions, particularly near 
100 K. (b) A magnetic order parameter (intensity of (1 -2 1) 
peak) that is influenced by moments aligned with c reveals 
a suppressed moment at low T. (d) Projection of the mag- 
netic moment per Fe onto the a6-plane and along c demon- 
strates the canting of the moment at 115 and 135 K. (e) The 
temperature dependence of the (10 5) peak reveals a struc- 
tural anomaly in the temperature region 100-140 K, which 
is primarily associated with corrugation of the Fe layers, 
(f) Anisotropic magnetic susceptibility^^ for single crystalline 
TlFei.6Se2 are consistent with a canting of the moments away 
from c between Ti ^100 and T2 ^UOK. 

Temperature dependent diffraction was utilized to 
study the origin of phase transitions at Ti ^ 140 K and 
T2 ~ 100 K, which were observed by measurements of 
the magnetic susceptibility (reproduced in Fig. |3]F), 
electrical resistivity, and heat capacity.!!^ Previously, 
we determined the c lattice parameter variation for a 
single-crystal using a 2-circle X-ray diffractometer, which 



showed behavior similar to that presented in Fig. 
Here, we show the whole temperature range for both lat- 
tice parameters determined with a four-circle neutron 
diffractometer, and the refinement results show a clear 
structural transition at ~100K. As observed in Fig. [3^,c, 
c increases abruptly near 100 K, while a decreases rapidly 
near 100 K and has unusual temperature dependence be- 
tween 100 and 200 K. No similar sharp anomaly was ob- 
served in unit cell volume at 100 K. The anomalous tem- 
perature dependence of a and c above ~200 K is likely 
due to the nearby order-disorder transition at 460 K. 

In the magnetically ordered state, the Fe layer is no 
longer flat but develops a slight corrugation with the 
block of 4 up spins (down spins) moving up (down). The 
parameter S is used to quantify this corrugation, and 
is given by the distance along c that separates the Fe 
spin blocks. In the absence of corrugation, Fel would 
be found at z=^. Therefore, 5/c is equal to twice the 
deviation from this ideal position of Fel. At 340 and 
250 K, refinements yield S ~ 0.04 A . The temperature 
dependence of the corrugation obtained by refining full 
diffraction data sets is shown in Fig. [4j The corruga- 
tion is only present in the magnetic state, as observed 
by a non-zero value below T/v, which indicates a strong 
coupling of magnetism to the lattice. The corrugation 
begins to increase rapidly near Ti, and reaches a maxi- 
mum value of S ^ 0.13 A at 115 K before a rapid de- 
crease at T2 to a value of S ^ 0.07 A , which is similar 
to that observed near room temperature. As observed in 
Table [T| the structural changes across Ti and T2 are dom- 
inated by changes in the lattice constants and changes in 
the corrugation. The position of Fel in Table [T] clearly 
demonstrates the corrugation, and the corresponding in- 
fluence on Sel (located above or below a spin block) is 
also be observed. 

The temperature dependence of the corrugation can 
be tracked via the intensity of the (10 5) peak, which 
partly reflects the corrugation and vanishes in the limit 
of zero corrugation. The (10 5) peak also gains intensity 
through a canting of the moment, which is discussed be- 
low. X-ray diffraction is therefore necessary to examine 
the extent to which the temperature dependence of the 
(10 5) peak is related to changes in the nuclear structure. 
Single crystal X-ray diffraction was performed at 200 K 
and 135 K, and reciprocal space maps with [1 0] cut are 
shown in Fig. |4j In these reciprocal space maps, one can 
clearly see the Bragg peak (10 5) and its equivalent (1 
-5) at 135 K but not at 200 K. To quantify the increase 
in intensity of (1 5) related to changes in the nuclear 
structure between 200 and 135 K, the intensity of (1 5) 
is integrated and normalized to the integrated intensity 
of (2 0), which is independent of temperature over this 
range (see Fig. The relative integrated intensities 

int(l 5)/int(2 0) are 0.43(8) at 200 K and 1.17(8) at 
135 K. This proves the nuclear contribution at (1 5) 
rises significantly between 200 K and 135 K, and thus the 
magnetic contribution to the data presented in Fig. |3^ is 
relatively small. 
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ferromagnetic system, the magnetic order parameter in- 
creases upon cooling below T/v (Fig. [sJd). Similarly, the 
magnetic susceptibility (x in Fig. [sjF) is relatively inde- 
pendent of T for fields {H) perpendicular to the spin 
orientation, while x decreases below T/v for H\\c . This 
behavior is maintained down to about 150 K, where the 
order parameter reaches a maximum and begins to de- 
crease. Interestingly, this corresponds to an increase in 
X for H\\c and a decrease in x for if±c . Upon cooling 
below ~100K, the rate of decrease in the order parame- 
ter increases, and x again behaves in a manner consistent 
with the observations above ~150K. The data for x are 
taken from Ref. 16 and were collected using a field of 5 T 
(zero-field cooled). 

At the temperatures above 150 K and below 100 K, 
the block- checkerboard antiferromagnetic model with the 
moments along c provides a good description of the neu- 
tron diffraction data. At temperatures between 100 K 
and 150 K, however, in-plane components of the moment 
are required. A canting of the moment is consistent with 
the trends in magnetic susceptibility shown in Fig. [3]F. 
The resulting block checkerboard and canted spin struc- 
tures are shown in Fig. |5j where only the Fe atoms and 
moments are shown for clarity. 



FIG. 4. (color online) Single crystal X-ray diffraction patterns 
with [1 0] cut showing the (10 5) and (1 0-5) Bragg peaks 
have much larger intensity at 135 K compared to 200 K. This 
is associated with an increased corrugation {5) of the Fe layer, 
which is obtained from refinements of the neutron diffraction 
data. 5 is defined as the distance along c between the spin 
blocks shown in Fig. l2jl. 



The phase transitions at Ti and T2 have clear struc- 
tural components. The transition at 100 K appears par- 
ticularly sharp, as evidenced by the jump in c (Fig. S]^) 
and the abrupt decrease in the corrugation (Fig. [4|, while 
the structural effects that begin at ~150K occur more 
gradually (Fig. |3|i,c). Both phase transition tempera- 
tures are well defined in the magnetic susceptibility data 
(Fig. [3]F), and the behavior of the magnetic structure 
through Ti and T2 are now considered. 

The (1 -2 1) peak shown in Fig. is utilized as a 
magnetic order parameter. This refiection is purely mag- 
netic, as indicated by its vanishing intensity at T^v, and 
is sensitive to the c component of the moment. The in- 
tensity of (1 -2 1) tracks that of the moment obtained by 
complete refinement, which is shown in Fig. [3]i. For the 
majority of temperatures, the total moment is equal to 
the moment along the c -axis (Mc). Upon cooling, this 
moment reaches a maximum at ~150K and then begins 
to decrease. A clear feature is observed near T2, below 
which the intensity of (1 -2 1) decreases rapidly. 

The magnetism is naturally broken into three tem- 
perature regions. Upon cooling below Tn ~430K, a 
block-checkerboard antiferromagnetic structure sets in 
with spins aligned along c . As expected for an anti- 




FIG. 5. (color online) (a) Block-checkerboard antiferromag- 
netic order with collinear Fe moments along c and (b) canted 
Fe moments as determined by neutron diffraction at 115 K. 

Through refinements of complete neutron data sets, 
the total moment is found to saturate to 2.3(1) /i^ at 
150 K and maintains this value until roughly 115 K, 
while the c -axis component of the Fe moment drops 
to 2.0(1) /iB • Between 100 K and 150 K, an a6-plane 
component of the moment {Mah) is observed. As shown 
in Fig. [3]i, Mah is only non-zero in the region between 
100 K and 150 K, which corresponds with the changes in 
X and the corrugation. The angle between the moment 
and the c -axis is 20(2)° at 140 K and 27(3)° at 115 K. 
Upon further cooling to 100 K, the block-checkerboard 
structure is recovered, but with a suppressed moment 
of 1.4(1) /i^, as shown in Fig. 3]i. Given the satura- 
tion of the vacancy ordering at high temperatures, this 
temperature-dependence of the moment is unlikely to be 
caused by changes in the degree of vacancy ordering. 

The magnetic ordering may be restricted to the va- 
cancy ordered regions, as proposed for the alkali metal 
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compounds However, the measurements presented 
here do not ahow the magnetic behavior of the vacancy 
ordered and disordered regions to be determined inde- 
pendently. 

Based on the results presented above, the tempera- 
ture dependence of the crystal and magnetic structure 
of TlFei.6Se2 can be summarized as follows. From 430- 
150 K the magnetic moments are ordered in the block- 
checkerboard arrangement with moments lying along the 
c direction. The Fe layer shows a slight corrugation in 
this temperature range. Upon cooling, at Ti ^ 140 K the 
moments begin to cant away from the c -axis. As the mo- 
ments cant, the corrugation of the Fe layer is significantly 
increased, indicating strong magnetoelastic interactions. 
This continues down to T2 ~ 100 K, at which point the 
canting abruptly vanishes and the higher temperature 
magnetic structure is recovered. However, the magnitude 
of the moment below T2 is significantly reduced from the 
value above Ti. These changes in the magnetism upon 
cooling through T2 are accompanied by two notable crys- 
tallographic changes: the corrugation sharply decreases 
and the c -lattice constant sharply increases, both reach- 
ing values similar to those seen above Ti. 

The observed increase in c and decrease in the mag- 
netic moment from a maximum of 2.3(1) /iB near 150 K to 
1.4(1) jiiB below T2 are remarkable. It is particularly un- 
usual to find a smaller magnetic moment at lower temper- 
ature with no difference in magnetic structure. Orbital 
ordering is predicted to occur in Fe layers of TlFei.6Se2 
when the vacancies are fully ordered, and provides a 
simple explanation for the observed behavior. In the 
absence of orbital ordering, each Fe has about two d- 
electrons distributed among three nearly degenerate t2g 



levels (with two lower lying, filled, Cg levels) giving a mo- 
ment of about 2 jj^B per Fe. In the orbitally-ordered state 
proposed by Lv et al.,'^the t2g orbitals are split so that 
each Fe has either dxz or dyz as the lowest energy level. 
Thus, orbital-ordering would decrease the occupancy of 
the dxy orbital as the electrons fill the lower energy states. 
The resulting increase in population of the out of plane 
{dxz and dyz) orbitals is consistent with the observed in- 
crease in c due to electrostatic repulsion. In addition, 
double occupation of the lowest lying t2g orbital would 
result in a decrease in the magnetic moment. With two 
t2g electrons per iron and complete orbital-ordering, the 
moment would be expected to vanish. However, the ob- 
servation of a non-zero moment at low temperatures is 
not surprising, since incomplete vacancy ordering would 
likely preclude perfect orbital ordering. 

This structural investigation of TlFei.6Se2 reveals that 
two low temperature phase transitions originate from 
competition between magnetic ground states and the 
coupling of magnetism to the lattice. The observation 
of competing ground states in TlFei.6Se2 highlights the 
complexity of the intercalated iron-selenide compounds. 
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